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We present cosmological perturbations of kinetic components based on relativistic Boltzmann equa- 
tions in the context of generalized gravity theories. Our general theory considers an arbitrary 
number of scalar fields generally coupled with the gravity, an arbitrary number of mutually inter- 
acting hydrodynamic fiuids, and components described by the relativistic Boltzmann equations like 
massive/massless collisionless particles and the photon with the accompanying polarizations. We 
also include direct interactions among fluids and fields. The background FLRW model includes 
the general spatial curvature and the cosmological constant. We consider three different types of 
perturbations, and all the scalar-type perturbation equations are arranged in a gauge-ready form 
so that one can implement easily the convenient gauge conditions depending on the situation. In 
the numerical calculation of the Boltzmann equations we have implemented four different gauge 
conditions in a gauge-ready manner where two of them are new. By comparing solutions solved 
separately in different gauge conditions we can naturally check the numerical accuracy. 
PACS numbers: 98.80.Cq, 98.80.Hw, 98.70.Vc, 04.62-Fv 



I. INTRODUCTION 

The relativistic cosmological perturbation plays a fun- 
damental role in the modern theory of large-scale cosmic 
structure formation based on the gravitational instabil- 
ity. Due to the extremely low level anisotropics of the 
cosmic microwave background (CMB) radiation, the cos- 
mological dynamics of the structures in the large-scale 
and in the early universe are generally believed to be 
operated as small deviations from the homogeneous and 
isotropic background world model. The relativisitic cos- 
mological perturbation analysis works as the basic frame- 
work in handling such cosmological structure formation 
processes. Recent observations of the CMB anisotropics 
in the small- angular scale by Boomerang and Maxima- 1 
experiments [1II4I, for example, confirm dramatically the 
validity of the basic assumptions used in the cosmologi- 
cal perturbation theory, i.e., the linearity of the relevant 
cosmic structures. 

Soon after the discovery of the CMB by Penzias 
and Wilson in 1965 §, Sachs and Wolfe in 1967 |] 
pointed out that the CMB should show the temperature 
anisotropy caused by photons traveling in the perturbed 
metric which is associated with the large-scale structure 
formation processes based on the gravitational instabil- 
ity. The detailed dynamics at last scattering is not im- 
portant in the large angular scale which can be handled 
using the null geodesic equations. Whereas, the physical 
processes of last scattering including the recombination 
process are important in the small angular scale where we 
need to solve the Boltzmann equations for the photon dis- 
tribution function |^ . When we handle the evolutions of 
the collisionless particles, like the massive/massless neu- 



trinos or the collisionless dark matters, we need the cor- 
responding Boltzmann equations as well. 

The relativistic gravity theory, including Einstein's 
general theory of relativity as a case, is a non-Abelian 
gauge theory of a special type. The original perturbation 
analysis was made by Lifshitz in 1946 based on Einstein 
gravity with the hydrodynamic fluid |^ . In handling the 
gauge degrees of freedom arising in the perturbation anal- 
yses in the relativistic gravity, Lifshitz started by choos- 
ing the synchronous gauge condition and properly sorted 
out the remaining gauge degrees of freedom incompletely 
fixed by his gauge condition. Other approaches based 
on other (more suitable) gauge conditions were taken 
by Harrison using the zero-shear gauge in 1967 Q and 
by Nariai using the comoving gauge in 1969 Each 
of these two gauge conditions completely removes the 
gauge degrees of freedom. Now, we know that the zero- 
shear gauge is suitable for handling the gravitational po- 
tential perturbation and the velocity perturbation, and 
the comoving gauge is suitable for handling the density 
perturbation. Since each of these two gauge conditions 
completely fixes the gauge transformation properties, all 
the variables in the gauge condition are the same as the 
gauge-invariant ones: that is, each variable uniquely cor- 
responds to a gauge-invariant combination of the variable 
concerned and the variable used in the gauge condition. 

The gauge-invariant combinations were explicitly in- 
troduced by Bardeen in 1980 |^; see also Lukash 1980 
]lOt for a parallelly important contribution. This be- 
came a seminal work due to a timely introduction of the 
early inflation scenario [ pH which provides a casual mech- 
anism for explaining the generation and evolution of the 
observed large-scale cosmic structures. We believe, how- 
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ever, a practically more important suggestion concerning 
the gauge issue was made by Bardeen in 1988 and 
it was elaborated in [|3|. In the gauge theory it is well 
known that proper choice of the gauge condition is of- 
ten necessary for proper handling of the problem. Either 
by fixing certain gauge conditions or by choosing cer- 
tain gauge- invariant combinations in the early calculation 
stage we are likely to lose possible advantages available in 
the other gauge conditions. According to Bardeen "the 
moral is that one should work in the gauge that is math- 
ematically most convenient for the problem at hand" . In 
order to use the various gauge conditions as advantages 
in handling cosmological perturbations we have proposed 
a gauge-ready method which allows the flexible use of 
the various fundamental gauge conditions. In this paper 
we will elaborate further the gauge-ready approach for 
more general situations of the generalized gravity theo- 
ries including components described by the relativistic 
Boltzmann equations. 

Our formulation is made based on the gauge-ready ap- 
proach; using this approach our new formulation of the 
cosmological perturbation is more flexible and adaptable 
in practical applications compared with previous works. 
Also, the formulation is made for a Lagrangian which 
is very general, thus includes most of the practically in- 
teresting generalized versions of gravity theories consid- 
ered in the literature. We pay particular attention to 
make the contribution of the kinetic components in the 
context of the generalized gravity theories. As an ap- 
plication of our gauge-ready approach made in this pa- 
per, we implemented the numerical integration of Boltz- 
mann equations for CMB anisotropies in four different 
gauge conditions. In addition to the previously used syn- 
chronous gauge (without the gauge mode) and the zero- 
shear gauge, we also implemented the uniform-expansion 
gauge and the uniform-curvature gauge in a gauge-ready 
manner. These two gauge conditions have not been em- 
ployed in the study of the CMB power spectra previ- 
ously. We will show that by comparing solutions solved 
separately in different gauge conditions we can naturally 
check the numerical accuracy. 

In §^ we present the classical formulation of the cos- 
mological perturbations of fields and fluids in the con- 
text of generalized gravity in a unified manner; i.e., di- 
verse gravity theories are handled in a unified form. The 
formulation is based on the gau ge-ready strategy which 
is explained thoroughly in j|ll E|. In §[II we present the 



gauge-ready formulation of the kinetic components based 
on the relativistic Boltzmann equations in the context of 
generalized gravity again in unified manner; i.e., we han- 
dle the massive/massless collisionless particles and the 
photon with Thomson scattering simultaneously, and all 
three-types of perturbations are handled in a single set 
of equations. In §IV we extend the formulation to in- 



clude the photon with polarizations, and implement the 
numerical calculation of the CMB temperature and po- 
larization anisotropy power spectra. Our present code 
is based on Einstein gravity including the baryon, cold 



dark matter (CDM), photon (including polarizations), 
massless/massive neutrinos, the cosmological constant, 
and the background curvature, both for the scalar- and 
tensor-type perturbations. The scalar-type perturbation 
is implemened using several gauge conditions where some 
of them are new. We explain how to generalize easily 
the Boltzmann code in the context of the generalized 
gravity theories including recently popular time varying 
cosmological constant. §^ is a discussion. In the Appen- 
dices ^ and ^ we present the conformal transformation 
properties of our generalized gravity theories and the ef- 
fective fluid quantities. In the Appendix ^ we present 
useful kinematic quantities appearing in the 3 -f 1 ADM 
(Arnowitt-Deser-Misner) formulation and the 1 + 3 co- 
variant formulation of the cosmological perturbation the- 
ory. 

We set c = 1. 



II. CLASSICAL FORMULATION 

A. Generalized gravity theories 

We consider a gravity with an arbitrary number of 
scalar fields generally coupled with the gravity, and with 
an arbitrary number of mutually interacting imperfect 
fluids as well as the kinetic components. As the La- 
grangian we consider 



y(</)^) 



(1) 



R is the scalar curvature. 0^ is the I-th component of 
N scalar fields. The capital indices I,J,K,... — 1, 2, 3, 
. . . , N indicate the scalar fields, and the summation con- 
vention is used for repeated indices, g = det{gab) where 
a,b, . . . are spacetime indices. f{(t>^, R) is a general alge- 
braic function of R and the scalar fields (j)^ , and gij{4>^) 
and V{(j)'^) are general algebraic functions of the scalar 
fields; f{(j>^,R) and V{(p^) indicate f {(t>\ . . . , (j)'^ , R) 
and V{(f>^ , . . . , (f)'^). We include a nonlinear sigma type 
kinetic term where the kinetic matrix gjj is considered 
as a Riemannian metric on the manifold with the coordi- 
nates (j)^ . The matter part Lagrangian Lm includes the 
fluids, the kinetic components, and the interaction with 
the fields, as well. 

Equation (|l]) contains many interesting gravity the- 
ories with scalar fields as subsets. Einstein gravity is 
a case of the minimal coupling with the gravity, thus 
/ = i?/(87rG); this case still includes the nonlinear sigma 
type couplings among fields, and for the minimally cou- 
pled scalar fields we have gu = Sjj. The general cou- 
plings of the scalar fields with gravity and the nonlin- 
ear sigma type kinetic term generically appear in various 
attempts to unify the gravity with other fundamental 
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forces, like the Kaluza-Klein, the supergravity, the super- 
string, and the M-theory programs; these terms also ap- 
pear naturally in the quantization processes of the grav- 
ity theory in a way toward the quantum gravity. The 
Lagrangian in eq. (^) includes the following generalized 
gravity theories as subsets [for simplicity, we consider one 
scalar field with (j) = (j)^ and gij — gii(0)]: 

(a) Einstein theory: / = -^^R, gii = 1, 

(b) Brans-Dicke theory: / = ^(l)R, 9ii = -g^, = 0, 

(c) Low-energy string theory: 

/ = e-^R, = -e-^ F = 0, 

(d) Nonminimally coupled scalar field: 

/ = (8^-?'^')^' 311 = 1' 

(e) Induced gravity: 



f^ec^-'R, gii = l, y-iA(^2_„2)2^ 
(f) i?2 gravity: / = ^{R + R^/QAP), 0-0, 



(2) 



etc. These gravity theories without additional fields and 
matters can be considered as the second-order theories. 
However, even with a single scalar field, the R) grav- 
ity is generally a fourth-order theory. Although such 
gravity theories do not have an immediate interest in the 
context of currently considered generalized gravity theo- 
ries, one simple example is a case with / = /i(0)/2(-R) 
where f2{R) is a nonlinear function of R. 

By the conformal transformation eq. ([^) can be trans- 
formed to Einstein gravity with nonlinear sigma model 
type scalar fields, and the transformed theory also be- 
longs to the type in eq. (|l|); see the Appendix Au- 
thors of 10 considered a less general form of Lagrangian 
than in eq. (|l|) in perturbation analyses; however, since 
they used the conformal transformation, they actually 
considered Einstein gravity with nonlinear sigma type 
couplings. 

Variations with respect to gab and (f>^ lead to the grav- 
itational field equation and the equations of motion: 



Gab ~ p 



Tal 



+ 5/J (^(^^a^b - 



{f^RF~2V)gab + F^a,b-9abF' 



8nGT^f\ 



(3) 



b'-'\ + -if-2vr+r 



rpb 

^a-b 



L 



where F = df/dR; 

^7 — 1 „IL 
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I. J 



91 J, ^JK 



(9 



LJ,K 



(4) 
(5) 

is the inverse metric of 
gLK.j - 9,jk,l), and Vj = 



dV/{d(t)'). Equation (|) follows from eqs. and the 

Bianchi identity. Tab is the energy-momentum tensor of 
the matter part defined as 5{^/^Lm) = ^y/—gT°'^Sgab- 
We have assumed the matter part Lagrangian Lm also de- 
pends on the scalar fields as Lm — Lm (matter, gaf,, 0^). 



In eq. (jij) term considers the phenomenological cou- 
plings among the scalar fields and the matter. In eq. 
(0) we introduced an effective energy-momentum tensor 

Tab^'' where the matter Tab includes the fluids and the 
kinetic components. The effective fluid quantities to the 
perturbed order are presented in the Appendix H. Using 



^ab 



we can derive the fundamental cosmological equa- 
tions in the generalized gravity without much algebra: we 
use the same equations derived in Einstein gravity with 
the fluid energy-momentum tensor and reinterprete the 
fluid quantities as the effective ones [|l5|. Direct deriva- 
tion is also straightforward. 

The matter energy-momentum tensor can be decom- 
posed covariantly into the fluid quantities using a normal- 
ized {u°'Ua = —1) four-vector Ua which is not necessarily 
the flow four- vector nm : 



Tab = fJ-UaUb 
h 



Phab + qaUb + qbUa + TTab, 
tl = TabU''u', p=^Tabh''\ qa=~T,dU^hi, 
T^ab = Tcdhahf - phab, 



(6) 



where hab = gab + UaUb is a projection tensor of the Ua 
vector, and qau'^ = = iTabU^, T^ab = T^ba, and tt" = 0. 
The matter energy-momentum tensor can be decomposed 
into the sum of the individual one as 



Tab^Y.^' 



{l)ab, 



and the energy-momentum conservation gives 



rp b 

^{i)a;b 



I 



Q 



{l)a 



(7) 



(8) 



where (i) indicates the z-th component of n matters with 
i,j,k,... = 1, 2, 3, ... n. The matters include not only the 
general imperfect fluids, but also the contributions from 
multiple components of the coUisionless particles and the 
photon described by the corresponding distribution func- 
tions and the Boltzmann equations. These kinetic com- 
ponents will be considered in §111 and IV. Q{^i)a takes 
into account of possible interactions among matters and 
fields. 



B. Perturbed world model 

We consider the most general perturbations in the 
FLRW (Friedmann-Lemaitre-Robertson- Walker) world 
model. As the metric we take 



2A) dTf ~ 2a^Badr]dx°' 



-a\g^^l+2C^0)dx''dx'', 



(9) 



where a{t) is the cosmic scale factor and dt = adrj. 
A{x,t), Ba{yi,t), and Ca/3(x, t) are generally spacetime 
dependent perturbed order variables. Ba, Ca, and Cap 
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are based on g^^, i.e., indices are raised and lowered with 

(3) 
9af3- 

The scalar fields are decomposed into the background 
and perturbed parts as 



(10) 



and similarly for R and F. In the following, unless nec- 
essary, we neglect the overbars which indicate the back- 
ground order quantities. 

The energy-momentum tensor is decomposed as 



a 



[qa + (M + P)Ua] = (a* + P)Va, 



(11) 



where Va and tt^'^^^ are based on g^^^. Va is a frame- 
independent definition of the velocity (or the flux related) 
variable . In the multi-component fluid situation from 
eq. (^ we have 



(12) 



and similarly for p, Sp, [p + p)va, and tt*'^''^ 



degrees of freedom for the tensor-type perturbation in- 
dicate the graviational wave. Whereas, two out of four 
degrees of freedoms, each for the scalar-type and vector- 
type perturbations, are affected by coordinate transfor- 
mations which connect the physical perturbed spacetime 
with the fictitious background spacetime. This is often 
called the gauge effect and a way of using it as advan- 
tages in handling problems will be described in § [I E . It 



is convenient to introduce the following combinations of 
the metric variables: 



X = a (/3 07) , I- 



3 {Ha - ifi) 



A 



a 



(14) 



where an overdot indicates a time derivative based on 
t, and H = d/a; A is a comoving three-space Lapla- 

cian, i.e.. Ax = x'"q. Later we will see that these com- 
binations are spatially gauge- invariant. The perturbed 
metric variables have clear meaning based on the kine- 
matic quantities of the normal-frame four-vector, see eqs. 

We introduce three-space harmonic functions depend- 
ing on the perturbation types. The harmonic functions 
based on g^^^ are introduced in PJl?]]: 



y(s)l7 



1 



C. Decompositions 

In the spatially homogeneous and isotropic background 
we can decompose the perturbed variables into three dif- 
ferent types, and to the linear order different perturba- 
tion types decouple from each other and evolve indepen- 
dently. We decompose the metric perturbation variables 
A, Ba, and Cap as 



A = a, 



(13) 



where (s), (v) and (t) indicate the scalar-, vector- and 
tensor-type perturbations, respectively. -Bi"\ and 
^q*/3 ^'"^ based on 5^^^, and a vertical bar indicates a 
covariant derivative based on g^^, for t(a/3) symbol see 
below eq. (p2|). The perturbed order variables a(x, t), 
P{yi,t), (p{x,t), and 7(x, i) are the scalar-type metric 
perturbations. Ba\'x.,t) and Ca\'x.,t) are transverse 
(^^"'"|q — — C^''^^°'\a) vector-type perturbations corre- 
sponding to the rotational perturbation. C^*2(x, i) is a 



transverse-tracefree (C^ a — — C ^|^) tensor-type 
perturbation corresponding to the gravitational wave. 
Thus, we have four degrees of freedom for the scalar-type, 
four degrees of freedom for the vector-type, and two de- 
grees of freedom for the tensor- type perturbations. Two 



-'q/3 — 



1 



is) 



Y 



1 



(3)y(s 



^al3 7 — ^ap^ 



^ap ^ 

vit) - V 

^ap — 



1. 



(v) 



(t) 
Pa 



-A - 

— ^ — ^ap ' 

(15) 



where k is a wave vector in Fourier space with fc = |k| ; the 
wave vector for individual type of perturbation is defined 
by the Helmholtz equations in eq. (15). In terms of the 
harmonic functions we have a(x, 77) = a(k, ry)!^'**^ (k; x) 

and similarly for /3, 7 and ip; B^"'^ = ci"' ee 



(v) 



and C. 



(t) _ 

ap 



c'-*-'y^^\ Since we are consider- 



ing the linear perturbations the same forms of equations 
will be valid in the configuration and the Fourier spaces. 
Thus, without causing any confusion, we often ignore 
distinguishing the Fourier space from the configuration 
space by an additional subindex. Also, since each Fourier 
mode evolves independently to the linear order, with- 
out causing any confusion we ignore the summation over 
eigenf unctions indicating the Fourier expansion. 

The perturbed scalar fields 6(j)^ in eq. (p^ ) only couple 
with the scalar-type perturbations, and are expanded as 



(5</.^(x,t) = 5?i^(k,t)r(")(k;x), 



(16) 



and are similarly for 6R and SF as well. 

Now, we consider perturbations in the fluid quantities. 



We decompose Va and tt 
bations as 



(3)a 



into three-types of pertur- 
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^(3)/3 ^ + + (17) 

The energy-momentum tensor in eq. ( pi] ) becomes 
T0 = -(/i + 5M), 

r° = {^l + p) ^^^^ + (m + p) 

= (p + ^p) + 7r(^)r(^>; + TT^^'^Y^^^; + 

(18) 

In terms of the individual matter's fluid quantities we 
have 



where fi, p and follow eqs. (|19y2 



Equations 



(19) 



and similarly for p, Sp, (/i -|-p)d'^'"\ and Tr^*''"'*^ We use 
the notation introduced by Bardeen in 1988 ||l2|; com- 
parison with Bardeen's 1980 notation [|j can be found in 
§2.2 of [l3| ; compared with our previous notation in fl^ ] 
we have TT^") = /a^)a and = -fc/[a(^ + p)]*. Wc 
often write v = w*-*^. 

The interaction terms among fluids introduced in eq. 
(^ are decomposed as: 

Q(i)o = -a + -4) + 5Q(i)\ , 

Q(,)„^jg¥(^) + jgVW. (20) 

From eq. (||) we have 

„ _ . -.1 
2^ Q{i) = r/(^ , 

I 

Y^ji^}=-TjS4>^, E4)=0. (21) 

Thus, RHS of the second equation in eq. (H) contributes 
only to the scalar-type perturbation. 

D. Background equations 

The equations for the background are: 
1 



3F 



M + IgiJ^'^' -^{f-RF-2V)- 3HF 



K 
7^' 



H = yn + p + gij 



F- HF 



K 



(22) 
(23) 
(24) 

(25) 
(26) 



^3|) follow from Gq and G?_— 3Gq components of Eq. 
|)^ respectively. Equation (^Sj) follows from eq. fl). 
Equation (|2^) follows from eq. (|^). By adding eq. ( pq ) 
over components we have 




/i + 3i?(^ + p) = T/t; 



(27) 



By setting F — 1/ {8ttG) we can recover 87rG factor in 
Einstein gravity. The gravity theory in eq. (^ includes 
the cosmological constant, A. The cosmological constant 
introduced in eq. (0) as an additional —A term can be 
simulated using either the scalar field or the fluid. Us- 
ing the scalar field we let V" —> V + A/{8ttG). Using 
the fluid, since A contributes = —Agat/ {8t^G) to the 
energy-momentum tensor, we let fi fi + A/(87rG) and 
p ^ p — A/ (SttG). This causes a change only in eq. (|2^). 
In the presence of the kinetic components we addition- 
ally have the Boltzmann equations for the components 
and the sum over fluid quantities should incl ude th e con- 
tributions from the kinetic components, sec §IIIC. 



E. Gauge strategy 

In the following we explain briefly our gauge-ready 
strategy. Due to the general covariance of the relativis- 
tic gravity theory we need to take care of the fictitious 
degrees of freedom arising in the relativistic perturba- 
tion analysis. This freedom appears because the rela- 
tivistic gravity is a constrained system: there exist some 
constraint equations with only algebraic relations among 
variables. In the perturbation analysis this is known as 
the gauge degree of freedom. The gauge freedom in the 
perturbation analysis arises from different ways of defin- 
ing the correspondence between the perturbed spacetime 
and the fictitious background. For example, by intro- 
ducing a spacetime dependent coordinate transforma- 
tion, even the FLRW background can be changed into 
a perturbed form which is simply due to the coordinate 
(gauge) transformation. Only in a special coordinate sys- 
tem the FLRW metric looks simple as in eq. (|^) without 
perturbations. 

Similarly as in other gauge theories, there are some re- 
dundant degrees of freedom in the equations which can be 
fixed without affecting the physics. Certainly it would be 
advisable, and is often essential, to take a proper gauge 
condition which either simplifies the mathematical anal- 
yses or allows an easier physical interpretation. Usually 
we do not know the best gauge condition (which differs 
depending on each problem) a priori, but it is desirable 
(actually often necessary) to find out the best one. In 
this regards, the advantage of managing the equations in 
a gauge-ready form was suggested by Bardeen in 1988 
]l2[ , and the formulation was elaborated in |13|. 

Contrary to many works in the literature which often 
refers the gauge freedom causing problems in the theory, 
we believe that, as in the other gauge theories (e.g., the 
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Maxwell theory and the Yang-Mills theory), the gauge 
freedom can/should be used as an advantage over solv- 
ing each specific problem. Our set of gauge-ready form 
arrangement of the equations, by allowing simple adop- 
tation of different gauge conditions, will allow the opti- 
mal use of the advantageous aspect of the gauge degrees 
of freedom present in the theory. To that purpose all 
the scalar-type perturbation equations are presented in a 
uniquely significant (see below) spatially gauge-invariant 
form but without fixing the temporal gauge condition. In 
this way, we can easily implement the several available 
temporal gauge conditions depending on the situation, 
and in this sense, the set of equations is in a gauge-ready 
form. The tensor-type perturbation describing the gravi- 
tational wave is gauge-invariant, and the vector- type per- 
turbation describing the rotation is presented using the 
uniquely significant gauge-invariant combinations of the 
variables. The particular choice of a gauge implies no 
loss of generality. If a solution of a variable is known in a 
specific gauge, the rest of the variables, even in the other 
gauges, can be easily recovered. Therefore, if possible, it 
would be convenient to start from the gauge condition 
which allows an easier manipulation of the equations. 
However, since the optimum gauge condition is usually 
unknown a priori, often it is convenient to carry out the 
analyses in the available pool of various gauge conditions 
and to find out the distinguished gauge condition; such 
analyses in the single component situations were carried 
out in the fluid ||l^, in the scalar field |jl^, and in the 
generalized gravity theories |^^. Our experience tells 
that different gauge conditions fit different problems, or 
even different aspects of a given system. Often, problem- 
atic aspects of the gauge freedom appear if one sticks to a 
particular gauge condition from the beginning and if that 
gauge condition turns out to be not a suitable choice for 
the problem. Our gauge-ready strategy is not a partic- 
ularly new suggestion in the contex of the gauge theory 
except that such a strategy, and its systematic use, has 
been largely ignored in the cosmology literature despite 
its rather apparent advantage. In the present work we 
extend the formulation in |13| to more general situation 
including kinetic components and arrange the equations 
for the convenient usage in diverse situations. 

In the perturbation analyses we have to deal with two 
metric systems, one is the physical perturbed model and 
the other is the fictitious background model. The gauge 
degrees of freedom arise because we have different ways 
of corresponding the perturbed spacetime points with 
the arbitrary background spacetime points. Since we 
are considering the spatially homogeneous and isotropic 
background the spatial correspondences (spatial gauge 
transformation) can be handled trivially: according to 
Bardeen Q "Since the background 3-space is homo- 
geneous and isotropic, the perturbation in all physical 
quantities must in fact be gauge invariant under purely 
spatial gauge transformations." We will show that only 
the variables [3, 7, b^'"\ and c'-"-' depend on the spatial 
gauge transformation. But these appear always in the 



combinations x a-nd ^P^") in eq. (|TJ) which are spatially 
gauge- invariant combinations; see eq. (|l|). These com- 
binations are unique in the sense that other combina- 
tions fail to fix the spatial gauge degrees of freedom com- 
pletely. Thus, using these (uniquely significant) spatially 
gauge-invariant combinations we take care of the effects 
of spatial gauge transformation of the scalar- and vector- 
type perturbations completely; the corresponding spatial 
gauge transformation properties of the kinetic compo- 
nents will be considered below eq. ([72|). 

Gauge transformation properties of the perturbed cos- 
mological spacetime were nicely discussed in 
pl|,p 17 1^,|2^. Under the gauge transformation of the 
form x°- — x'^ + the metric and the energy-momentum 
tensor transform as 



(28) 



thus, 



gabix") = gabix") - gab,cC - 9bcCa " 9acC,bi (29) 

and similarly for Tj^. By introducing ^' = (0 ~ rj) 

and = ^,a+^a\ with based on g^^j and ^i"^'" = 0, 
the perturbed metric quantities and the collective fluid 
quantities change as: 

a^a-i\ ^ = ip-H£_\ /3 = /? - 1^* -f a^, 

a 

v = v- -C, 
a 

Bi^) ^ B^:^ + ai^:\ C^:) -i^-), (30) 

and v''^\ ^fp-i Tr'*'"'*^ are gauge invariant. Thus, from 
eq. ( |l^ we have 



(31) 



and these are spatially gauge-invariant. From the scalar 
nature of (f)^ , R, F, and Tj we have: 

s^i = - , -^f , = sTi - Tie , 

SR = SR-R^\ SF = 6F-F^K (32) 

From eq. ( j^dj ) we notice that the tensor-type pertur- 
bation variables are gauge- invariant. For the vector- type 
perturbation we notice that ^^"^ defined in eq. (|lj) is a 
unique gauge-invariant combination. Thus, using ^("^ 
the vector- type perturbation becomes gauge-invariant. 
For the scalar-type perturbation using x instead of /3 
and 7 individually, all the variables are spatially gauge 
invariant. Considering the temporal gauge transforma- 
tion properties we notice several fundamental gauge con- 
ditions based on the metric and the energy-momentum 
tensor: 
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O 1 

Synchronous gauge: 


a = 


0, 


Comoving gauge: 


V 1 k 


= 0, 


Zero-shear gauge: 


X = 


0, 


Uniform-curvature gauge: 


if> = 


0, 


Uniform-expansion gauge: 


K = 


0, 


Uniform-density gauge: 


X . , - 
Ofl = 


— n 


Uniform-pressure gauge: 


Sp E 


eO, 


Uniform-field ((^^) gauge: 




= 0, 


Uniform- i? gauge: 




= 0, 


Uniform- _F gauge: 


SFi 


= 0, 



(33) 

etc. The names of the gauge conditions using x, ip and k 
can be justified: these variables correspond to the shear, 
the three-space curvature and the perturbed expansion 
of the n ormal-frame vector field, respectively, see eqs. 
(|C3|,|CT^). v/k = is a frame-invariant definition of the 
comoving gauge condition based on the collective veloc- 
ity. 

The original definition of the synchronous gauge in |^] 
fixed /3 = as the spatial gauge condition in addition to 
a = as the temporal gauge condition. In such a case, 
from eq. ( ^0| ) we notice that the spatial gauge fixing also 
leaves remaining (spatial) gauge degree of freedom. By 
using the spatially gauge-invariant combinations x ^-nd 
V we can avoid such an unnecessary complication caused 
by the spatial gauge transformation which is trivial due 
to the homogeneity of the FLRW background |l2) . From 
eq. ([l^) X is the same as af3 in the 7 = gauge condition. 
But in the P — gauge condition we have x = fl^7i thus 
7 is undetermined up to a constant (in time only) factor 
which is the (spatially varying) remaining gauge mode. 

By examining eqs. ( p0| , ^ , p2p^ ) we notice that, out 
of the several gauge conditions in eq. (|3^), except for 
the synchronous gauge condition, each of the gauge con- 
ditions fixes the temporal gauge mode completely; the 
synchronous gauge, a = 0, leaves spatially varying non- 
vanishing ^* (x) which is the remaining gauge mode even 
after the gauge fixing. Thus, a variable in such a gauge 
condition uniquely corresponds to a gauge-invariant com- 
bination which combines the variable concerned and the 
variable used in the gauge condition. Several interesting 
gauge-invariant combinations are the following: 



5fj,y = Sfj,- -fiv, (p^ = ip-Hx, 
^v = V- —V, 54)^ = S(j) - — 



k 

V X, 



LP \ 



H 



VS4>'- (34) 



For example, the gauge-invariant combination (50^ is 
equivalent to dcj)^ in the uniform-curvature gauge which 
takes ip = Q as the gauge condition, etc. In this way, 
we can systematically construct various gauge-invariant 
combinations for a given variable. Since we can make sev- 
eral gauge-invariant combinations even for a given vari- 
able, this way of writing the gauge-invariant combination 
will turn out to be convenient practically. 



In the multi-component case of fluids there exist some 
additional (temporal) gauge conditions available. From 
the tensorial property of T(^i)ab and using eq. ( p9| ) we can 
show 



(35) 



and vl^l .(^) 
tional temporal gauge conditions are 



TT/f^''''*-' are gauge invariant. Thus, the addi- 



Sn(^i)=0, (5p(,,) = 0, vi^i)/k = 0, 



0, 



(36) 



etc. Any one of these gauge conditions also fixes the 
temporal gauge condition completely. From the vector 
nature of Q(^i)a and using eq. (|^) we have 



SQ{i) 



r ' — T 



(v) 



is) 



J, 



is) 



(37) 



As mentioned previously, in general we do not know the 
suitable gauge condition a priori. The proposal made in 
Ip^ , ^ is that we write the set of equation without fixing 
the (temporal) gauge condition and arrange the equa- 
tion so that we can implement easily various fundamen- 
tal gauge conditions. We call this approach a gauge-ready 
method. Any one of the fundamental gauge conditions in 
eqs. (|3^ , |36|) and suitable linear combinations of them 
can turn out to be a useful gauge condition depending 
on the problem. A particular gauge condition is suitable 
for handling a particular aspect of the individual prob- 
lem. The gauge transformation properties of the kinetic 
components will be co nside red in § [II ; see paragraphs 
surrounding eqs. (|7^,|98,105). 



F. Scalar-type perturbation 

In this section we present a complete set of equations 
describing the scalar-type perturbation without fixing 
the temporal gauge condition, i.e., in the gauge-ready 
form. 

Definition of k: 



1 1 fc2 
tp = Ha- + Q^X- 
3 3 



(38) 



ADM energy constraint (Gg component of the field equa- 
tion): 



k^-3K 







hi 









SfJ. + gij<t> 



<?W0''-(/-2U)_^ 
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-3HSF + [3H + 3H^ -'^ ] SF 



(39) 



Energy conservation of the fluid components [from 
^(j)0;6 = Q(*)o and using eq. dH)]: 



Momentum constraint (G" component): 

3 [a 



-3K 3F 



{^1+ p)v + gijcj) 



+6F~HSF . 
ADM propagation (G^ — ^S'^G^ component): 

X+(^H+^^X-a-^^ji^^n^^^+SF 

Raychaudhuri equation (G^ — Gq component): 



K+ 2iJ H kH d 

\ 2F 2F 



(40) 



(41) 



3H 



2F \ 



6F + 3HF + Agij, 



= ^U^i + 3Sp + Agij 



+ 



2gij,K^'4>' + if - 2V) 



K 



+3SF + 3H5F + \ -6H^ + '^—^^ ] SF 



Scalar fields equations of motion: 
{2V-fy\^+V'j^^^, 



^ (k + d) + (20^ + 3H4,^ + 2r5^ 



+ ^F-'SR-Sr^. 
Trace equation (GJJ component): 

SF + 3H5F +(^-^\5F+ ^gu^'S^^ 



guM'^' + 2 if- 2V)^ 
{Sfi - 35p) + F{K + d) 



1 



K 



^gij^'<p-' + 2F + 3Hf) a - ^F6R. 



Scalar curvature: 
6R = 2 



k - AHk + { ^ - 3H ] a 



k^~3K 
-2— 



(42) 



(43) 



(44) 



(45) 



+ +P(i)) K + (5(9(,). 



(46) 



Momentum conservation of the fluid components (from 



1 



+P{i) 



. 2k^-3K (,) ^ 

3 — p — '^w 



(47) 



By adding properly eqs. (ji^:^) over all components of 
the fluids, and using properties in eqs. (19, 2l| , |2^ ), we get 
equations for the collective fluid quantities as: 



Sfi + 3H (^/i + 8p) = (/i + p) K - 3Ha v 

\ a 



(48) 



1 



a\^^+p) 

_ k 
a 



[a'^{fi+p)v] 



1 



fj.+p 



(49) 



It is convenient to introduce 



Spik, t) EE cl{t)S^iik, t) + e(fc, t), S=^, 



;(t) 



P 



(50) 



Equations (^8[]49|) provide a redundantly complete set 
for handling the most general scalar-type perturbation of 
the FLRW world model allowed by the Lagrangian in cq. 
(|l|); for example, eq. (||) follows from eqs. (||,||,|5l). 
Equations (^,|9|) follow from eqs. (|^,|^. Following the 
prescriptions below eq. ( [27| ) these equations also include 
the cosmological constant; an introduction of A term does 
not appear explicitly in our set of equations in the form 
eqs. (H-H). Notice that eqs. (^|9|), which follow 
from fluid energy- momentum conservation in eqs. (^J^), 
are not affected formally by the generalized nature of 
the gravity we are considering. In §[II and ^ we will 
see that the presence of kinetic components additionally 
introduces the corresponding Boltzmann equations, and 
their contributions to the energy-momentum content can 
be included as the individual fluid quantity in the above 
set of equations. 

Equations (|3^p9|) are written in a gauge-ready form. 
In handling the actual problem we have a right to impose 
one temporal gauge condition according to the mathe- 
matical or physical conveniences we can achieve. As long 
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as we choose a gauge condition which fixes the temporal 
gauge mode completely, the resulting equations and the 
solutions are completely free from the gauge degree of 
freedoms and variables are equivalently gauge-invariant. 
Some recommended fundamental gauge conditions are 
summarized in eqs. ( ^3|j3^ ). Equations (|3^^ are de- 
signed so that we can easily accomodate any of these 
gauge conditions. 

If we take an ansatz for term in eq. (11) as 



to the perturbed order we have 
-J 



(51) 



(52) 



where we used w° = ^(1 — a). Such a phcnomcnological 
damping term was considered in 



in the Boltzmann equations though, see eqs. (p5,102). 
Equations (^,^ which are independent of the field 
equations tell that in a medium without anisotropic stress 
terms 7r|j^^ and the mutual interaction terms among com- 
ponents J(i) , the angular momentum combination of in- 
dividual component is conserved as 



Angular Momentum ■ 



X a X V 



(i) 



constant in time. 



(57) 



The presence of anisotropic pressure can work as the sink 
or the source of the rotational perturbation of the in- 
dividual fluid. Conservation of the angular momentum 
combination of the rotational perturbation in Einstein 
gravity was noticed in the original work by Lifshitz |6j . 

In the presence of kinetic components we additionally 
have the corresponding Boltzmann equations, and the 
components contribute to the anisotropic pressure in the 
above equations, see §g and 0. 



G. Rotation 

The equations for the vector-type (rotational) pertur- 
bation are: 



fc2 - 2K 



1 r 



a'' {fj, +p) 



2K 



2k 



2K 



2k 



(v) 



(53) 
(54) 

) 

(55) 



Equations (53) follows from G° component of Eq. (||), 
and eq. (Ei) follows from T^.j, = 0. Equation ( |55| ) follows 
from eq. (p). By adding eq. ( |55| ) over all components we 
have eq. ^J). Notice that eqs. ( [5^j55| ) are not affected 
formally by the generalized nature of gravity theory. In 
fact, these two equations are derived from the conser- 
vations of the energy- momentum tensors in eqs. (|5|J^) 
without using the gravitational field equation. The pres- 
ence of kinetic components additionally introduces the 
corresponding Boltzmann equations, and contributes to 
the fluid quantities in the above equations, see §111 and 

The vorticity tensors based on the frame-invariant 
four- vectors are (see the Appendix ^) : 



(i)a/3 



(56) 



Thus, we have w = \fuj°\jabj2 and similarly for '^(i)- 
Equations ( |53| - |5^ ) show that the fluid velocities of the 
rotational perturbation do not explicitly depend on the 
generalized nature of the gravity, whereas, only the met- 
ric connected with the rotation mode, \E''^''\ depends 
on the nature of generalized gravity; term appears 



H. Gravitational wave 

The tensor-type perturbation (gravitational wave) 
equation in Einstein gravity was derived originally by 
Lifshitz in We can derive easily the wave equation 
for the most general situation covered by the Lagrangian 
in eq. (|^) as 




2^ (t) _ 1 Vtt^*) 
/ 



(58) 

which follows from component of eq. (|^) using eqs. 
(|l^,^,^ . The generalized nature of the gravity appears 
in the F terms: one in the damping term and the other in 
modulating the amplitude of the fluid source term. This 
equation is valid for the general theory in eq. (^), and the 
presence of arbitrary number of the minimally coupled 
scalar fields (with general g/j) does not formally affect 
the equation for the cosmological gravitational wave. The 
presence of kinetic components additionally introduces 
the corresponding Boltzmann equations, and contributes 
to the anisotropic pressure in the above equations, see 
Q and |l^. 

Equation (|58|) can be arranged in following form 



2K 



(t) 

(0' 



Vt 



= aVFc'-*\ zt = aVF, 



(59) 



where a prime denotes the time derivative based on rj. 
In the large-scale limit, thus ignoring the k^ term in eq. 
(p9|), and assuming = and tt*^*) = 0, we have a general 
integral form solution ll^ 
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c'-*\k,t) =c{k)~d{k) 



-dt. 



(60) 



where c(fc) and d{k) are integration constants for rela- 
tively growing and decaying solutions, respectively. This 
solution is valid considering the general time evolution 
of the background dynamics as long as the perturbation 
is in the superhorizon. The growing solution is simply 
conserved in the superhorizon scale and the generalized 
nature of the gravity does not affect the conserved nature 
of the growing solution. Only in the decaying solution the 
generalized gravity nature appears explicitly. 

Similar equations and solutions as above can be derived 
for a single component scalar-type perturbation in unified 
forms for the fluid, the field, and the generalized gravity 
theory as well |23|. 



III. KINETIC THEORY FORMULATION 

A. Relativistic Boltzmann equation 

The evolutions of coUisionless particles and the photon 
are described by specifying distribution functions which 
are governed by the corresponding Boltzmann equations. 
The relativistic Boltzmann equation is given as p3| 



dX 



f 



dx" df 

'dXdx^ 

C[fl 



dp'' df 
'dXdp^ 



(61) 



where /(a;°,p^) is a distribution function with the phase 
space variables a;" and = dx°'/dX, and C[f] is the col- 
lision term. The energy-momentum tensor of the kinetic 
component with mass m is given as 



T('^4 = / 20(/)<5(p=p, + m2)pV/^/^dV 



0123 



(62) 



Assuming the mass-shell condition, after integrating over 
p^, we have 



t; 



ab 



^d^p 



'3^123 



bol 



(63) 



Equations (^||j|) together with eqs. (61,6|), including 
T(c)ab in the individual fluid energy-momentum tensor, 
provide a complete set of equations considering the con- 
tribution of a component based on the distribution func- 
tion (we call it a kinetic component). The corresponding 
fluid quantities can be identified using eq. (^). In the 
case of multiple kinetic components, we have eqs. (0^) 
now valid for the individual kinetic component. The cor- 
responding fluid quantities of the individual component 
can be identified using eqs. (^,0). 



B. Boltzmann equation in perturbed FLRW 

Under the perturbed FLRW metric in eq. (||), using 
p°- as the phase space variable, eq. ( |6l|) becomes 



+ 2 f Aa - 



a 



2-pp +T'^J^ pPp-fj — ^df]. 



dpo 
(64) 



In handling the Boltzmann equation and the energy- 
momentum tensor in perturbed FLRW spacetime, it is 
convenient to introduce a special phase space variables 
based on a tetrad frame. In literature we find several 
different choices for the phase space variables p5|- p7| , p7t ■ 
As the phase space variables we use {q, 7") introduced as 



„o 



p" = — (1 — A)^/q'^ + rn^a^, 

-,(3) 



(65) 



where 7" is based on with 7"7q = 1. The advantage 
of this choice in our gauge-ready approach will become 
clear below eq. (^). Using (9,7") as the phase space 
variables eq. (|64| ) becomes 



2^ y ^ ^7" ) 



\/ q^ + m^a^ 



\/ q^ + m?a- 



-:{l + A)C[f]. 



q— 
dq 

(66) 



We decompose the distribution function into the back- 
ground and the perturbed order as 



fin, a:", g, 7") = /(r?, q) + 5f{ii, x", g, 7"). 



(67) 



Assuming that the collision term has no role to the back- 
ground order, in which the Thomson scattering is the 
case, we have 



/' = 0. 



(68) 



Thus, / is a function of q only. The energy-momentum 
tensor in eq. (p3) becomes 



(69) 



rj^ab [ a bf q^dqdVLq 



^q^ + rn^a^ 
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The fluid quantities defined in eq. (|11|) become: 
q^dqdflq 



P(c) 



3a4 



(A^Cc) +P(c))v(c)a 



a* 



Sf-faq^dqdVtq, 



'(c)/3 



a' 



Under the gauge transformation = x"" + ^""j consid- 
ering p'^ = dx'^/dX, we have = + jp**. Using the 
definition of q in eq. (pa), and using eq. (pOj) we have 



(71) 



Since / depends only on q, we have (5/ = 5f — {df/dq){q- 
g), thus 



dq \ q a ' 



(72) 



Notice that with our phase space variables in eq. ( |65| ) 
the perturbed distribution function Sf is spatially gauge- 
invariant. Thus, our choice of the phase space variables is 
particularly convenient for the gauge-ready formulation 
where, as a strategy for later convenient use, we do not 
fix the temporal gauge condition while fixing the spa- 
tial gauge condition without losing any advantage; our 
Sf is spatially gauge invariant. We can show that the 
gauge transformation property of Sf is consistent with 
the gauge transformation properties of the fluid quanti- 
ties identified in eq. (|70|). 



C. Background equations 



Equations (|22|-|27D describe the evolution of the FLRW 
world model. The sum over fluid quantities in eq. ( |l9| ) 
should include the kinetic components. To the back- 
ground order, from eq. ( [70| ) we have 

fJ-ic) = ^ I ftq^dq, = ^ I f'^d.q, (73) 



where e{q,'ri) = \/ q^ + m?a?; hereafter, the mass m ap- 
pears only in e. Thus, for massless particles we have 



-J / /9 dq, = -H[c 



(74) 



We can show that eq. ( pq ) applies to the kinetic com- 
ponents as well with (z) — (c) both for the massless and 
massive particles. This identification gives 



Qic) - 0. 



(75) 



In the case we have the matter (m), radiation (r), and 
massive neutrino {vm), it is convenient to introduce 



17a 



Mc 

A 
ZIP' 



Mc 



8^' 



K 
(76) 



The matter includes the baryon and the CDM with = 
f2f, -|- and the radiation includes photons and the 
massless collisionless particles like massless neutrino {v) 
with Q,r = + r^i/. In such a case we have 



M = X! '^(O ^ + + Pl-I + f^V + l-i-Vrr^ + MA, 

; 

P^^P(l) "^P-t+Py+Py^ 
/ 



(77) 



where we have recovered the cosmological constant using 
the prescription below eq. (p7|). 

The distribution function of the fermi/bose (± sign) 
particle is given as 



/(e) 



1 



(78) 



where is the number of spin degrees of freedom, and 
hp and fc^ are the Planck and the Boltzmann constants. 
If the decoupling of the massive particle occurs while it is 
relativistic^hus for neutrino mass much less than IMeV, 
e in eq. (78) can be approximated as g, and afterward 
the distribution function is well approximated by Fermi- 
Dirac distribution with zero rest mass. 

Since = 0, from eq. (26) we have /im = Mb + 
Mc Of a^^ and /i^ = M7 + M^ °f Foi^' the photon and 
massless neutrinos we have 

1 1 7/4\'/' 
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1/3 



(80) 



where is the number of massless neutrino species. 
Equations (^-26) together with eqs. ([75|,[77|) describe 



the background evolution in the context of generalized 
gravity theories. The generalization to include multi- 
component massive/massless collisionless components is 
trivial: we simply consider the fluid quantities in eqs. 
(F1,F3) for each component. 



D. Massless particle 

For a massless particle, m = 0, it is convenient to 
introduce a frequency integrated perturbed intensity 



SI{x'',^) = AQ{x'',^) 



^ JSfq'dq 
jfqHq ■ 



(81) 
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For photons, unless we have an energy injection process 
into the CMB, the spectral distortion vanishes to the 
linear order. Assuming photon's distribution function 



/ 



9s 



1 



(82) 



we can expand the temperature fluctuation 8 = 5T/T in 
the following form as well 



Sf 



.91 
'dq 



(83) 



In terms of 6 the perturbed part of eq. (pq) becomes 



d 



d 



+collision term. 



The fluid quantities in eq. ([70|) become 



M(c) 



J3)a 



(85) 



Using the spatial and momentum harmonic function 
introduced by Hu et al in p8[ , we expand 



e(x,,y,7)^EE E e[J(k,^)G[y(k;x,7) 



k i=0 m=-2 



(86) 



GjJ(k;x,7)^(-*) 



'sjl^Yr{Kl)e''^''-^\ (87) 



where i > |m|, and to = 0, ±1, ±2 correspond to the 
scalar-, vector-, and tensor-type perturbations, respec- 
tively. Thus, for the scalar-type perturbation we have 



G 



(0) 



(-j)^P,(k.7)e'^^(''^-^), 



where (5(k; x) is a spatially dependent phase factor which 
depends on the harmonic functions, see eq. (p2|); in the 
flat background we have e*^^'''''' = e*'''^. As the normal- 
ization we have 



^(m)^(m')* 



1 



47r 2£ + 1 
We have the recursion relation 
d 



(89) 



d 



(m) 



2£- 



where k! 







1 and for ^ > 1 



'(^2„m2) 



1- (£2 



1 



H) 



(91) 



The harmonic functions are introduced such that we 
have (0), (±1), and (±2) superscripts instead of (s), (w), 
and {t) indices introduced in eq. (15). In terms of the 
spatial harmonic functions we identify 



Mm) ai o 

^(|m|)-7 ••■7 



Qi .. .ai 



(92) 



Then, from eqs. (^,|T^ we can show: 



G 



(0) 
(0) 



' (1) ~ ' a ' 



2 /l_3i£G'") 



(84) g(±i) ^^ay(±i) ^./i„2— G^*'^ -T^^'^r^*'^- 
g[^^^ ^^'^^PY^f\ (93) 



Now, from eqs. ( |84| ) using the expansion in eqs. 
and the recursion relation in eq. (pG) we can derive 



(m) 



(94) 



where M^^^ and G^^^* are the metric perturbation and 
the collision term, respectively. The collision terms for 
the Thomson scattering in photon distribution function 



The metric perturbations in eq. (^4|) can be calculated 
using eqs. (f3^,P3h as 



together with polarizations will be considered in §IVA 

'i 

/ 
V 



7I\/i-2ff.i.(±i) 







a(±2) 



(95) 



where the rows indicate to = 0, ±1, ±2, and the columns 
indicate I — 0, 1,2, respectively. Using eq. ( p6|) the per- 
turbed order fluid quantities in eq. ( p5[ ) become: 

^l^(c) (0) 1 



M(c) 



_e(o) 



(1)' 



(c) 
M(c) 



= 9 



(0) 
(2)' 



V 



(±1) _ q(±1) 



TT 



(±2) 



(1) 



M(c) 



15 



=9 



1 - 2 



(±1) 

(2) ' 



— 9(±2) 

M(c) " 15 ''(2) ' 



(96) 
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where we used eqs. ([1 



Under the gauge transformation, from eq. (72) we have 

e = e + He + -ec.i"- (97) 

a ' 

By expanding ^* = Ek^*^'°' 

we have 



e 



(0) _ 0(0) 



(0) - ^(0) ' ' ^(1) - ^(1) Q 

:>("') 



(98) 



and other are gauge-invariant. These are consistent 

with the identifications in eq . (|9q) and the gauge trans- 
formation properties in eqs. ( |3C|p5[) . Thus, the temporal 

gauge conditions fixing o||j| and o|'j'| can be considered 
as belonging to eq. (|3^). 

E. Massive collisionless particle 

For massive collisionless particles the perturbed Boltz- 
mann equation in eq. (pq) becomes 



oq 



(99) 



In the massive case we can expand directly the perturbed 
distribution function. Instead of 6f we use the following 
variable 



From eqs. (70 100 ) the perturbed order fluid quantities 
become: 

(M(c) +Pic))v^,^ = 0(1) [-9^)1 dq, 

(to = 0,±1). 



,(±1) 



1-3 



\/3 Stt /■ ^(±1) f df \ q- 



^_2K lba*J (2) 1^ 



dq, 



(±2) Stt f ^(±2) / df\ q^ 



(103) 



where we used eqs. ( |l7yiqj93| ). In the massless limit, 
assuming is independent of g, the fluid quantities 



in eqs. (^,103) reduce to the ones in the massless case 

eqs. (^4|j96l) with e[y = ej;"^ 



Under the gauge transformation, using eq. ( |72| ) we 
have 



^ ^ e + He + --e,ai' 

q a 



(104) 



Thus, 



0(x,?7,(?,7) 



Sf 



^Y.T. E e[,";)(k,,y,g)G[™)(k;x,7). (100) 

k /=0 m=-2 

From eq. (p9|), using eqs. (100 90,1^,|9^) we can derive 



g k 



e 



^ ' e ( 



1 



21 



1 



2£ + 3 



rn A(m) 



M, 



(m) 



(101) 



7fVl-2ff*(^^) 
-c(±2) 



(102) 



where i > \m\. The rows and columns of M^^^^ indi- 
cate m = 0,±1,±2, and / = 0,1,2, respectively. In the 
massless limit M^^^ reduces to M^"^^ in eq. (p5|), and 



eq. ( 101) reduces to eq. (|94[). Thus, we can regard eqs. 
( |101| , 102) as being valid for both massless and massive 
collisionless particles. 



;(o) 



~(0) 



(0) 



(1) 



q a 



(105) 



and other 8 



(m) 



These are consistent 



/(.^-j are gauge-invariant, 
with the identifications in eq. (|103|) and the gauge trans- 
formation properties in eqs. (| 30| , p5D . Thus, the temporal 



and g|'j'| 



can be considered 



gauge conditions fixing 6|pj 

as belonging to eq. (|3^). 

Thus, we have complete sets of perturbation equations 
for three-types of perturbations including a single com- 
ponent massive collisionless particle. The eq. (101) to- 
gethcr with the gravitational field equations in eqs. (p8[- 



3H55|J58|), and the fluid quantities for the massive par- 



ticle in eq. (103) provide the complete sets. In the case of 
multi-compone nt massiv e collisionless particles we simply 
consider eqs. (101, 103| ) for each component of massive 
collisionless particle. The equations for the scalar-type 
perturbation are designed in a gauge-ready form. The 
collective (or the sum over individual) fluid quantities 
in eqs. (|39|- ^ , fi[|54 ,|5l ) include the kinetic components, 
whereas (i) in eqs. (46|^4^; jSj) does not include the kinetic 
components. Using eqs. ( 103|,[75|), h owever, we can show 
that eq. (IT0l|) gives eqs. MaBmsl) with (i) = (c). This 



identification gives 
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SQ{c) = = J| 



(s,v) 



(106) 



This follows because we have assumed coUisionless situa- 
tion (and with no direct interaction between the kinetic 
component and the other components). F or t he case of 
photon with Thomson scattering, see eq. (121). 



IV. CMB ANISOTROPY 
A. Thomson scattering and polarizations 

Besides the photon distribution function for the tem- 
perature (or total intensity) fluctuation / = /e, we have 
three other photon distribution functions describing the 
state of polarization, /g, fjj, and /y. 0, Q, U, and 
V form the four Stokes parameters. We will ignore the 
fourth Stokes parameter V describing a circular polar- 
ization because it cannot be generated through Thom- 
son scattering in standard FLRW cosmological models. 
While the temperature behaves as a scalar quantity Q 
and U do not. It is known that the combinations Q ±iU 
behave as the spin ±2 quantities |29|. Thus, while 
can be expanded in ordinary spherical harmonic Yim, 
Q ± iU should be expanded in the spin-weighted har- 
monics ±2Yim @|l|- 

Following the convention in p^ ] we expand & and Q ± 
ill in terms of the spin-weighted spatial and momentum 
harmonic functions: 

2 

e(x,r;,7)^5] ^ 5]ejJ(k,^)oG[™)(k;x,7), 

k m=-2 i 

Q(x,?7,7) ±if7(x,77,7) 

-EE EN?(k.'?)±<7'(k,^) 



k m=-2 



x±2G[^"^(k; x,7), 



(107) 



where 



(108) 



with oY;" = Y^"; thus oG["^ = G'-^^K We have a recur- 
sion relation |^ 

Mm) a _ { m Mm) _ m ^(m) \ 

(109) 



. inns (rn) 



where 



n= VP -I- (1 -f ImDX, 



(110) 



thus, compared with eq. ( |91| ) we have n x qk™ = /sk™. In 
the hyperbolic (negative curvature) background we have 
the supercurvature (0 < fc < ^|i^|) and subcurvature 
(fc > ^|iir|) scales for the scalar-type perturbation; by 
considering n > we exclude the supercurvature scale, 
H. It is convenient to have -^Y^™ = (-l)^^y^'", and 
other useful relations can be found in ||3^,|2^ . 
In terms of a notation 



T 




the Boltzmann equation can be written as 
7^ + -Tja7" = M[f]+C[f], 



(111) 



(112) 



where the metric perturbations and the collision terms 
are expanded as: 



m^T. E E^^, 

k m=-2 f 

^[0]-E E E^, 



(m)^(m) 
(t) ^(£) ' 



(m)^(m) 



(113) 



(114) 



k m=-2 



The collision and the polarization terms are not affected 
by the perturbed metric, thus M[Q ± ill] — 0. 

The collision term is derived using the total angular 
momentum method in eqs. (25,26) of |^ 



c[f] = -ff{n) 



J&f^+rV(b)a 
T I 





10 



(115) 



m=-2 



where is the baryon's perturbed velocity variable, 

and p(™) is given in eq. (52) of [Q; f = rieX^o-T where 
He is the electron density, Xe is the ionization fraction, 
and (Tt is the Thomson cross section. The time evolu- 
tion of determined by the recombination history. 
Using p(™) in |33] the collision term becomes 



C, 



(m) 



(m) 



^(0) 






p(0) 



-It)'' P'^'' 
P(±2) 



(116) 



C[Q±^U] = -rY: E E (47^ ± 

k m=-2 I 



(m) 



+V6P(")<5,2)±2G 



(m) 
(2) ' 



where 



p(m) ^ J_ (Mm) /pp, 
~ 10 V (2) ^Dii., 



(m) 
(2) 



(117) 



(118) 
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From eqs. (112,107), using the rec ursio n relation in eq. 
( |l09| ) and the collision term in eq. (117), we can show: 



(m) 
(^) 



n 

a 1^2^ 
2m 



1 



a £(^ + 1 



1 



E, 



(m) 



2^ + 3 



(119) 



(^) 



a 1^2^ 
n 2m 
"a £{£+!) 



(m) 
(l-l 

p(m) _ ■pirn) 



1 



2^ + 3 



(120) 



where £ > 2 and m > 0, and we have B 



(0) 



that with an identification B 



(-|™l) 



0. Notice 

(|m|) ^(±|m|) 



and oIm'™'"* satisfy identical equations for both signs 



'{i) 



Equations for e[™^ follow from eqs. (|94|,|95|,|116|). Polar- 
ization properties of the CMB in perturbed FLRW world 
model have been actively studied in the literature; some 
selected references are |3l|j3^ , |35| , p9| , ^ -|38|| . 

Now we have complete sets of equations for three- types 
of perturbations including Thomson sca ttered photons 
with polarizations. Equations (^ 95 llC) for the inten- 
sity (temperature) and eqs. ( 119 12C| ) for the photon 
polarization together with the gravitational field equa- 
tions in eqs. (|3^ ^ , ^3| - |55| , p8|) , and the fluid quantities 
for the massless particle in eq. ( p6[ ) provide the complete 
sets. For a massless collisionless particle eq. ( p^ re- 
mains valid with vanishing collision terms and the polar- 
izations. The generalization to include multi-component 
massless collisionless particles is trivial: we simply con- 
sider eqs. (^,^) for each component of massless colli- 
sionless particle. We can also include additional multi- 
comp qnent ma ssive collisionless particles by considering 
eqs. ( |l01 103) for each component of massive collision- 
less particle. The equations are designed in a gauge-ready 
form. 

The collective (or sum over individual) fluid quantities 
in eqs. 
Using eqs 
gives eqs. 



44 



54 5^) include the kinetic components. 



96,75), however, we can show that eq. ^ ^ 

(46, 47, HI) with (i) = (c); we are using the 
index c to indicate the kinetic component. This identifi- 
cation implies 



SQ{c) 
(±1) 



= 0, 



7(0) 4 a , 



J, 



/±1) 

'(c) 



(±r 



(121) 



Due to Thomson scattering there exists an interaction 
between photons and baryons. From eq. (UW we have 

E; ^Q{i) ^Q{b) + SQ(j) = and J^i 



J, 



(m) 
(7) 



0. Thus 



6Q(b) = 0, J^ 



(0) 

(fc) 



7(±1) 



= - J, 



(±1) 
(7) 



4 a I 
.,(±1) 



'(7) 



(±1) 
\b) 



- J\ 



ib) 



For the baryon (b) we have = = thus 
W(^b) — 0. However, we keep the sound speed of the 
baryon fluid which behaves as [BM 



2 _ P{b) _ _ 1 dlnTb \ ksTb 



M(b) 



(123) 



where u is the mean molecular weight. Thus, eqs. 
(1^,1^,11) become: 



Vib) 



■Jib) 



Hi 



(124) 
(125) 



For the cold dark matter (C) we additionally have Q(c) 
= ^Q(c) and J^^^j = 0, thus eqs. ( p6| , p7| ) become: 



3(C) 



k 

--V{C) 



3Ha + K, 



k 



Hvic) = -a- 
^ ' a 



(126) 
(127) 



We emphasize that compared with previous works, be- 
sides the equations are valid in the context of generalized 
gravity theories, our sets of equations in this paper are 
all in gauge-ready forms. 



B. Tight coupling era 

In the early universe, the Thomson scattering term is 
large enough and the baryons and the photons are tightly 
coupled. For large value of f(= = X~^) compared 
with H{= t^^ = X„^) and k/a{^ SttA^^), it is difficult 
to handle eqs. (^,125) numerically; the polarizations are 
negligible in that stage. In such a case, it is convenient 
to arrange the equations in the following way ]39| ] . From 
eq. (|125[) and the i — 1 component of eq. d94) we have 



l + r 



-Hv(^b) - 



^(b)<'ib) 



'\i) 



'(7) 



1 



Oo(0) 



5 '^20(2) 



(128) 



where we introduced r = |/i(.^)//i(f,). 
component of eq. (B4) we have 



From the 



e 



(0) 
(2) 



10 



e 



(0) 
(2) 



«2«(7) 



(3) 



4/1 - 3^1 i 



V6 (0) 
9 (2)' 



(129) 



where Cr = H/t. Thus, O^"! is of the c^-- or c^-^- 



(122) order higher compared with and we have Q f^l-i 
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^C t6|° |. From the £ = 1 component of eq. ( p^ ) and 
eq. ( |l28|) we have 



+ - ( cf,,6, 



a 



Cr 1 

Hl + r 
1 



1 



(fc)"(b) 



Oo(0) 



5 "-2^(2) 



(130) 



Taking a time derivative of eq. (|130|) and using eqs. 
(128,130), to the first order in expansion, we can derive 



1 



2r 
1 + r 



H{v 



(b) ~ ^(7) J 



Hl + r 
k 



(1 + ry 



(b) 



Ha - C(h^5, 



{b)0(b) 



-HS, 



(7) 



(131) 



where we assumed the radiation era with a cx t^/^, which 
apphes for Einstein gravity with neghgible K and A con- 
tributions in that era. Thus, Cr — Hcr, r = —Hr, and 



we used (c?^^)■ = —Hcff^^ which follows from eq. (122) 



-(b)' - " 

assuming Tb = T^. 

To the zeroth order in c^, from eqs. (12S.131) we have: 



2r 



'J{b) 



1 2r^ 

- 7TX7T2 ^ (^(7) - ^^(6)) 



1 + r 
k ( 1 



^- a + 

a V 1 + r 



(l + r)2 

4A^) + I + ,4 



(7) 



(132) 



Equations (125,131) imply V[b) ~ ^^(7) = to the zero-th 
order in c^. Thus, using v^b = ^^(b) = W(^), and ignoring 
the c^fc-)<^(6) term, we have 



v^b 



1 k ( r 1 

Hv^b H — [ a + -( 

1+r ^ a\ 1+r 4 



'(7) 



(133) 



Equation (124) and the i — component of eq. ( |94| ) give 

k 



\b) 



-Vjb — 3Ha + K, 



4 f k 

S(^) = — v~.b — iHa + K 

6 \ a 



(134) 
(135) 



and for the CDM we have eqs. ( |l26|Jl27l ). 

Thus, in the tight coupli ng era, instead of eqs. 
(^,125), we can use and eqs. (133- 135[). As the criteria 
for the tight coupling era we can use [14G|] 



Z> Ztc = 2000, Cr < 0.01, Cr < 0.01. (136) 

aH 

where 



X, 



14.47 / a 

ao 



H 
H~o' 



(137) 



If any one of the three conditions is violated we use the 
full set of equations based on the Boltzmann equation. 



C. Numerical implementation 

The different gauge conditions available in the scalar- 
type perturbation provide a useful check of the numerical 
accuracy |^ . The gauge-ready formulation is especially 
suitable for handling the case. When we solve the scalar- 
type perturbation equations we have a right to choose 
one temporal gauge condition. Any one of the funda- 
mental gauge conditions in eqs. ( |3^j3^ ) would be a fine 
gauge condition; except for the synchronous gauge con- 
dition any one of the other gauge conditions completely 
fixes the temporal gauge and the remaining variables are 
equivalent to the gauge-invariant ones. If we have the 
solution of a variable in a given gauge we can derive so- 
lutions of the rest of variables in the same gauge, and 
from these we can derive all the solutions in other gauge 
conditions. For such translations the set of equations in 
the gauge-ready form is convenient. Meanwhile, in the 
numerical study, if we solve a given problem in two dif- 
ferent gauge conditions independently, by comparing the 
value of any gauge-invariant variable evaluated in the two 
gauge conditions we can check the numerical accuracy. 

In the literature, the synchronous gauge is the most 
widely adopted gauge condition. The synchronous gauge 
does not fix the gauge condition completely. We can 
choose any of the gauge conditions mentioned in eqs. 
(33,36,38 105) as well. The comoving gauge condition 
closely resembles the synchronous gauge in the matter 
dominated era; however, see the various possible combi- 
nations of the comoving gauge conditions in eqs. ( |3^ ) 
available in the multi-component situation. The authors 
of adopted a comoving gauge condition which fixes 
the velocity variable based on the cold dark matter, thus 
V(^Q-^/k = 0. From eq. (127) we notice that — im- 
plies a — 0, the synchronous gauge condition; in the syn- 

" ^ which leads 



chronous gauge, however, we have V(^c) ^ 
to the remaining gauge mode. Thus, the synchronous 
gauge with an additional condition Ui-,^) = is equiva- 
lent to the ti((7) = gauge. We may emphasize that our 
equations in the gauge-ready form are ready to be imple- 
mented using any of these available gauge conditions. 

In the following, as an example, we consider Einstein 
gravity limit without fields but with arbitrary numbers 
of fluids and kinetic components; including the fields 
and the generalized gravity will affect the gravity sector 
only. In the numerical work we implemented the comov- 
ing gauge based on CDM (w(c) = which includes the 
synchronous gauge), the zero-shear gauge, the uniform- 
curvature gauge and the uniform-expansion gauge. The 
latter two gauge conditions were not previously used in 
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the hterature. These four gauge conditions fix the per- 
turbed metric variables. Under these gauge conditions 
we can see that the differential equations can be set up 
only using the variables which r epresent the fluids and ki- 
netic components: eqs. (101,34,119,120, Mprh provide a 



set of differential equations to be solved. The remaining 
metric variables can be expressed in terms of the fluids 
and kinetic quantities. The metric variables a, Hx 
(or y/a), and k/H (or qk) are dimensionless. Using eqs. 
( |38| - |42| ) we can express the metric variables in terms of 
the fluid quantities in each of these gauge conditions. 

In order to check the numerical accuracy we can eval- 
uate any gauge-invariant variable in two different gauge 
conditions; if the integration has good numerical accu- 
racy the gauge-invariant combination evaluated in all 
gauge conditions should give the same value. In order 
to have the same solution we need to start from the 
same initial condition. Thus, we need to have relations 
of variables among different gauge conditions. The set 
of equations in a gauge-ready form is convenient for this 
purpose. In the following we consider relations between 
the zero-shear gauge and the uniform-curvature gauge as 
an example. Using the gauge transformation properties 
in eqs. ( p0| , pl|j3^ ,|9g ,105) we can construct the following 
relations 

= ^(i) + 3(1 + w^,))ip = (5(j)^ + 3(1 + w^i))(p^, 
k 



aH 



fix: 



e(o) 



0(o!x + ^- 



e(o) 



and 9 



(0) 



(138) 



follows from eqs 



6|^j^ for £ > 2; and similarly for '■Px 



fix 



k^-3K 
-Hx^- 



E 



(139) 



Using eqs. ( 138| , 139D we can translate the solutions (in- 
cluding the initial conditions) in the zero-shear gauge into 
the ones in the uniform-curvature gauge, and vice versa. 

Each of the four gauge conditions considered above 
flxes the metric variables as the gauge condition, and 
uses the fluid, fleld and the kinetic variables as the un- 
known variables to be solved. In such cases we can make 
the numerical code which allows to choose one of the 
gauge condition as an option. In the numerical study it 
is known that the zero-shear gauge has numerical diffi- 
culty in setting up the initial condition at early universe 
1^,^. The other three gauge conditions show no such 
a difficulty encountered in the zero-shear gauge and run 
equally well. 

We have implemented our gauge-ready formulation 
into the numerical code. The code includes the baryon 
(&), CDM (C), photon (7), massless (i^) and massive neu- 
trino (vm) species, the spatial curvature, and the cosmo- 
logical constant. We included photon polarizations. We 



solved separately the gravitational wave with the accom- 
panying tensor-type photon intensity and polarizations, 
massless and massive neutrino species. The set of differ- 
ential equations is solved directly. For the recombination 
process we adopt the RECfast code by which is re- 
cently available in public. As the initial conditions, we 
have implemented the flve different non-decaying initial 
conditions available in the four component (&, C, 7 and 
v) system in the radiation dominated era p4| : these are 
the adiabatic mode, the baryon isocurvature mode, the 
CDM isocurvature mode, the neutrino isocurvature den- 
sity mode, and the neutrino isocurvature velocity mode, 
where the last one appears due to the kinetic nature 
of the neutrino perturbation. The complete set of ini- 
tial conditions was recently presented by Bucher et al in 
[Q in the = gauge condition. The correspond- 
ing initial conditions in the other gauge conditions can 
be obtai ned by u sing the gauge transformations similar 
to eqs. (138 139 ). In the early radiation dominated era 
we used the tight couplin g approximat ion for the baryon 
and the photon in eqs. (133 134,135) with the criteria 
in eq. (136). The code is complete in the context of 



Einstein gravity. We made no artiflcal truncations for 
multipoles of kinetic components in the photon intensity 
and polarization, and massless/massive neutrinos. For 
a useful truncation scheme, however, see Since the 
higher i multipoles are generated from the lower multi- 
poles, we monitered the values of highest multipoles of all 
the kinetic components and increased the allowed multi- 
poles automatically ||4^ . In this manner we included the 
quantities with higher multipole as long as the values are 
larger than certain minimum threshold value. We made 
the code so that we can choose the gauge condition out 
of the four different gauge conditions mentioned above as 
an option; we could try other gauge conditions as well. 
The inclusion of the scalar flelds and generalized grav- 
ity is a trivial generalization affecting only the gravity 
sector, and will be considered in future occasions. 
At the present epoch we have 



k 



2998 k 



aoHo 



ao 



-Mpc, 



(140) 



where Hq = lOOhkm/ secAIpc. Thus, k/ao = \Mpc~^ 
corresponds to Aq — 27rao/fc — 2TrMpc, and fc/(aoi?o) = 
2998//i; the comoving wave number k is dimensionless. 
In Einstein's gravity from eq. ( |2^ ) we have K/{aH)^ — 
il — 1. Thus, in the non-flat case we have 



ao 



2998Mpc 

^\i-no\h' 



(141) 



In the hyperbolic model, the curvature scale corresponds 
to /c^ = —K = 1 and the subcurvature scales (fc^ > 1) 
correspond to fc/ap > ^/l — floh / {2998Mpc) . In the 
spherical model, the wavenumber n introduced in eq. 
dllOD takes integers n = 3, 4, 5, . . . (n = 1, 2 corresponds 
to pure gauge modes [^). 
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In the following we present several results from our 
numerical study. In the numerical integration of the dif- 
ferential equations we adopted the Runge Kutta method. 
The integrations were made in the equal interval of In a. 
In Figures ^ (a,b) we showed the evolution of S{i)v^i-i 
which is the density perturbation of (i)-component in 
the corresponding comoving gauge condition based on 
the component 



= +3— (l + w(,))i;(,), 



(142) 



which follows from eqs. ( p5| , p6|) ; = P(i)/n(i) and 
we ignored The component (i) includes (6), (7), 

(C), {v) and (i^m)- The small scale considered in Fig- 
ure |l| (a) crosses the horizon in the radiation dominated 
era. After the perturbation comes inside the horizon, 
the baryon, photon, and massless neutrino show oscilla- 
tions, and after the recombination near log(a/ao) ~ —3 
the baryon decouples from the photon and catches up 
the evolution of the cold dark matter. The behavior of 
massive neutrino is also shown. The large scale pertur- 
bation in Figure ^ (b) crosses the horizon in the matter 
dominated era and the oscillations do not appear. 

In Figures |^ (a,b) we presented the behavior of ip^ as 
well, ipv was first introduced by Lukash in 1980 
see also Q, and is known to be one of the best con- 
served quantities in the single component situation: it 
is conserved independently of changing gravity theories 
or field potential in the super-horizon scale |p^ , po| , and 
independently of changing equation of state in the super- 
sound-horizon scale p8| ; for a summary, see |23|. In the 
Figure it shows nearly conserved behavior while in the 
super-horizon scale and in the matter dominated era after 
the recombination; however its amplitude changes near 
horizon crossing, and is affected by the recombination 
process if the scale is inside horizon. We showed detailed 
behaviors of (fy and ipy^.^ in Figures I and |. Figure | 
shows behaviors of if in various comoving gauge condi- 
tions based on fixing v or U(j-) for three chosen scales. We 
showed the evolution of ip^ for several different scales in 
Figure ||. 

In Figures || (a,b) we presented the evolution of 
^(i^m)v(^ ) which is the density perturbation of the mas- 
sive neutrino in the comoving gauge based on the massive 
neutrino. In Figure ^ (a) we considered a model domi- 
nated by the massive neutrino, showing the coUisionless 
damping of the neutrino density fluctautions . The 
case with substantial amount of the cold dark matter is 
presented in Figure ^ (b). In the massive neutrino the 
fluid quantities include the integral of the distribution 
function over the momentum variable, q in cq. (103): in 
our numerical work we considered about 100 values of q 
for a range of q/ (kBaoTo). 

As the wavenumber k increases, i.e., as we consider the 
smaller scale, we need to solve the larger number of the 
differential equations. As examples, for k/ao = 0.001 and 
O.lMpc^^ considered in Figures |l| (a,b) i is excited up 



to around 600 and 5000, respectively. We increased £ au- 
tomatically by monitoring the value of individual kinetic 
component (including the polarizations). 

Aspects of the roles of massless neutrino in the evolu- 
tion of cosmic structures were studied in |^ . The roles 
of massive coUisionless particles (massive neutrino is the 
prime example) as the hot dark matter in the context of 
structure evolution have been investigated in the litera- 
ture ||2|,|^,|5|M|,|9|^,||]. Further analytic studies can 
be found in p7|. The gravitational instability using the 
particle distribution function was originally studied by 
Gilbert in 1965 in the Newtonian context Ml. 



D. CMB anisotropy 

The anisotropics of the temperature can be derived by 
expanding the observed temperature in the sky into a 
spherical harmonic function as 



e(x,77o,7) = ^ af^(x)y^„(7). 



(143) 



The polarization anisotropics can be expanded in terms 
of the spin- weighted harmonic functions as 

Q(x,?7o,7) ±iC/(x,?7o,7) 

e 

E K™(x)±*aL(x)]±2r^™(7). (144) 

I m=-l 

We can derive: 



1 ^ 

Y E («^war™(x))x 

2 

(145) 



(2^-M) 



m=-2 



where X and Y can be any one of Q, E and B. In 
the flat and hyperbolic background {K < 0) we have 
n > 0, see below eq. ( |110| ). In the background with 
positive curvature, we have discrete n with n = 3,4, . . .; 
see below eq. (141). In such a case the integration should 
be changed to a sum over n with n = 3, 4, . . .. Due to the 
parity, we have Cf^ = = Cf ^ If the distributions 
are Gaussian, all statistical informations are contained 
in the three angular power spectra and one correlation 
power spectrum between Q and E: 



C, 



ee 



EE 



r<BB 



(146) 



Both the scalar-type and gravitational wave contribute to 
the correlation functions ®, Cf^ and Cf^, whereas 
only the gravitational wave (and the rotation) con- 
tributes to Isl. 
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In Figures I (a,b) and Fi gure 1^ we presented the power 
spectra of the scalar, and tensor-type perturbations. In 
both the scalar- and tensor-type pert urbation spectra, 
for the integration over k in eq. (|l45| ), we took 500 fc's 
in the equal interval of log A; for a range from fc/oo = 
10"-* - 0.5Mpc-^; to have Ce complete to ^ - 2000 
(which would well cover the MAP and Planck Surveyor 
results) we need kmax/o-o ^ 0.2Mpc^^. The spectra are 
filtered using a smoothing method. 

Pioneering work concerning the CMB anisotropy based 
on relativistic gravity and Boltzmann equation was made 
by Peebles and Yu in 1970 Early theoretical works 
can be found in p0|-|5^. Significant progress has been 
made in the theoretical side of the CMB anisotropy im- 
mediately following the first detection of the quadrupole 
and higher-order multipole anisotropics by the COBE- 
DMR and subsequent ground based experiments; some 
selected references are |3|^0[|3|,||,||,||,||,|||| . Au- 
thors of have developed a CMBfast code which cal- 
culates the CMB angular power spectra in an efficient 
way using the line of sight integration of the Boltzmann's 
equation. CMBfast is based on the synchronous gauge 
and is applicable to Einstein gravity. Authors of p2] 
have modified the code adopting a comoving gauge con- 
dition based on the velocity variable of the cold dark 
matter (our v^c)/k — gauge). In our code, we solve the 
full Boltzmann hierarchy without any approximation. A 
run takes less than an hour (without the massive neu- 
trino) in the workstation equally well for all three gauge 
conditions we used. 



The power spectra in eq. (146) are known to be sen- 
sitive to various combinations of the background world 
models (these include the Hubble constant, spatial cur- 
vature, cosmological constant, and density parameters of 
baryon, CDM, massless and massive neutrinos), the ini- 
tial amplitudes and spectra of both the primordial den- 
sity and gravitational wave, and the possible reionization 
history, etc. Thus, in return, the observational progress 
in determining the power spectra can give strong con- 
straints on the above mentioned parameters with higher 
precisions ||57[| . 

There has been a significant improvement of the CMB 
power spectrum measurements in the past decade, and 
further improvements are expected from the ground 
based, balloon, and flight experiments, and particularly 
from the planned MAP and Planck Surveyor satellite 
missions with a high-accuracy and small angular res- 
olution. The recent balloon observations of CMB by 
Boomerang and Maxima- 1 experiments have already pro- 
vided a strong constraint on the global curvature of our 
observed patch of the universe: the location of first peak 
in Figures ^ (a,b) corresponds, in models with f^o near 1, 
to -^peak — 200/-\/no whcrcas the Boomerang experiment 
shows £peak = 197 ±6, thus supporting a flat universe [0. 

The small I plateau region in Fig. || can be inter- 
preted as reflecting the primordial scale-invariant spec- 
trum, which has arisen from the quantum fluctuations 
in the context of the inflationary scenario. The small I 



corresponds to the large angular scale, and the plateau 
region corresponds to the super-horizon scale in the last 
scattering epoch where the local scattering would be 
unimportant. Thus, the Sachs- Wolfe effect based on the 
null-geodesic equations is expected to be enough to ex- 
plain the physics (based on the relativistic gravitation 
relating the spacetime metric with matter). Meanwhile, 
the oscillatory features in the large t (small angular scale) 
come from regions well inside the horizon at the last scat- 
tering, thus in addition to the gravity the local physics 
including direct couplings between photons and baryons 
is important. Now, the physics behind these oscillatory 
features is well understood as being due to the oscillatory 
evolution of the photon fluctuation (and tightly coupled 
baryons) pictured/frozen at the last scattering epoch: the 
oscillatory evolution of each k mode (reaching the last 
scattering epoch with different phase) together with the 
initial spectrum is reflected into the corresponding os- 
cillatory feature in k space which can be converted into 
the oscillatory feature of Ci in the I (angular) space. In 
hindsight, the original prediction of this oscillatory fea- 
ture can be traced back to Sakharov as early as 1966 
(which was before the discovery of CMB); for a clear 
exposition, see p9|-|6l[|, and for more elaborated forms 
see Q. More visually, the oscillatory feature in k space 
can be found in Peebles and Yu in 1970 [|| and Sunyaev 
and Zeldovich in 1970 and more developed analysis 
was made in Doroshkevich et al in 1978 |Q. Up to our 
knowledge, however, the first clear and complete picture 
of the oscillatory features (including the polarization as 
well as the isocurvature case) in Fig. |^ can be found in 
Figure 7 of Bond and Efstathiou (1987) H]. 



V. DISCUSSIONS 

Compared with the previous works we have made some 
notable advances in the formulation. The formulation 
is made for the general form of the Lagrangian in eq. 
(|l|) which is more general compared with our previous 
works. The kinetic theory treatments in §111 and |rv| are 
presented in a gauge-ready form for the scalar-type per- 
turbation. Also, the kinetic theory formulation is made 
in the full context of the generalized gravity theory cov- 
ered by the Lagrangian in eq. (0). 

For the scalar-type structure all the equations are ar- 
ranged in a gauge-ready form which enables the optimal 
use of various gauge conditions depending on the prob- 
lem. Usually we do not know the most suitable gauge 
condition a priori. In order to take the advantage of 
gauge choice in the most optimal way it is desirable to 
use the gauge-ready form equations presented in this pa- 
per. Our set of equations is arranged so that we can 
easily impose various fundamental gauge conditions in 
eqs. (33[3^), and their suitable combinations as well. 
Our notation for gauge-invariant combinations proposed 
in eq. (p4) is practically convenient for connecting solu- 
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tions in different gauge conditions as well as tracing the 
associated gauge conditions easily. 

In handling the Boltzmann equations numerically, 
we showed uniform-expansion gauge and the uniform- 
curvature gauge could also handle the numerical inte- 
gration successfully. By comparing solutions solved sep- 
arately in different gauge conditions we can naturally 
check the numerical accuracy. It may deserve examin- 
ing the physics of CMB temperature and polarization 
anisotropics in the persepective of these new gauge con- 
ditions and others which might still deserve closer look. 
Our set of equations in a gauge-ready form is particu- 
larly suitable for such investigations where we can easily 
switch our perspective based on one gauge condition to 
the other. 

In this paper one can find the general cosmological per- 
turbation equations which are ready for use in diverse 
FLRW world models based on the gravity theories in eq. 
(|l|). More attention will be paid on the generalized ver- 
sions of gravity theories especially in the context of early 
universe in future. In such a context, the formulation 
made in this paper would be useful for studying the struc- 
ture formation aspects of the future cosmological models. 
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APPENDIX A: CONFORMAL 
TRANSFORMATION 

By the conformal transformation the gravity theories 
included in eq. (0) can be transformed into Einstein 
gravity [ ]62| . Under the conformal transformation of the 
spacetimc metric, gab = ^^gab, and the field redefinition, 
n = VSttGF EE e5\/I'>, eq. (|l|) becomes 



with 



V = 



1 



{l6T:GFy 



i2V-f + RF) , 



(Al) 



(A2) 



where we have ignored Lm term. Thus, in general, since 
■0 = ■0(0^, i?), we have an additional minimally coupled 



scalar field -0. However, \i ip = ip{(j)^) which is the case 
for / — F{if>^)R and for the gravity theories in eq. (||), 
eq. (Al) becomes 



1 



167rG 



R 



where 



1 



91 J 



1 



SvrG \F 



-giJ + 4>,i4',J 



(A3) 



(A4) 



Relations we need to derive the above results can be 
found in (l|j6|. 

As a simpler situation we consider a case with gu = 



and — ip{4') where i,m = 2, 3, 
introducing 



, N, and = 



1 (9n 
8nG V F 



d(j)'^ + dip' 



we can show that eq. ( Al ) becomes 
1a 1 1 



L 



IGttG 



R 



SttGF^ 



By 
(A5) 

(A6) 



where we have a canonical form kinetic term of (p. Equa- 



eq. ( 



tion (AG) als o fo llows directly from eq. (AS). Notice that 
eqs. (Al,A^,A6) all belong to our original Lagrangian in 

(IT 



The conformal transformation in the context of cos- 
mological perturbation has been considered in ||lj,|l^,^ . 
We decompose the conformal factor il into the back- 
ground and the perturbed part as 



i7(x,<) = n{t) 

Thus, we have: 
n = VSttGF = e^^'f', 



2F 2 



(A7) 



-S,p. (A8) 



In [ [I5y63| ] we have shown that the only changes under the 
conformal transformation are the following 

a = aQ, di = ^dt, a = a + ^O, ip ^ ip + 5VL. (A9) 

Thus, in our multicomponen t si tuation, assuming the 
conditions used to derive eq. (|A6|) are met, we have: 




56 66 6F , 

^^i^-F- ^^^^^ 

J. t 



[In ||15|,|63[ we considered the situation with a single field 
with gii = Lo{4>). In the present case gn and gim are 
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arbitrary algebraic functions of (j) and cf/ ] From these 
we can also show that 



H 



(All) 



are invariant under the conformal transformation. Re- 
lations among (j), 0, and F in the individual gravity are 
summarized in Table 2 of . The advantages of using 
the conformal transformation in cosmological perturba- 
tion as a mathematical trick to simplify the analyses are 
presented in |15|,|63|. 



APPENDIX B: EFFECTIVE FLUID QUANTITIES 

We present the effective fluid quantities based on the 
effective energy-momentum tensor introduced in eq. (|^). 
The effective energy-momentum tensor in eq. (||) is de- 
composed into the effective fluid quantities similarly as 
in eqs. To the background order we have: 



(off) 



1 r 
F 



-3HF 



8nGp 



(off') _ 



1 



F 

F + 2HF 



if~RF~ 2V) 



\U-RF-2V) 



„(cff) _ n _ (off) 



(Bl) 



The scalar-type perturbed order effective fluid quantities 
are [use eqs. (B4,B5) in poi]: 



1 



-3HSF + [ 3H + 3H^ - -^j 6F 



-Fk+ (3HF- gij<j> 



-4 (giJ^^S^'' + SF- H5F - Fa 
a r \ 



F 



:{SF~Fx) 



(B2) 



The vector-type effective energy-momentum tensor is: 
F 



F 1 

'F2^ 



I" I ( 



I" 



(B3) 



The tensor-type effective energy-momentum tensor is: 
^(5T^'^l-FC$). (B4) 



APPENDIX C: KINEMATIC QUANTITIES 

The 3-1-1 ADM equations and the 1-1-3 covariant 
equations are convenient in analysing the cosmologi- 
cal perturbations, ||9|, |6^ , ^ , |65t . The kinematic quantities 
and the Weyl curvatures appearing in the formulations 
are useful to characterize the variables used in the per- 
turbation analyses. In the following we present various 
quantities appearing in the two formulations in the con- 
text of our perturbed FLRW metric. For the basic sets 
of the ADM and the covariant equations, see §VI in 
|l6| and the Appendix in . 

The covariant decomposition of the normalized 
{n°'na = —1) normal-frame vector field Ua provides clear 
meanings of the perturbed metric variables. The normal- 
frame vector field is introduced as 



no = -a(l + A), 



0. 



(CI) 



The kinematic quantities based on the normal-frame vec- 
tor are 



na;bn , 



(C2) 



8^G(5p(°ff) = -hp + gjji 



gij,K(b'r + {f~2V) 



S(t>- 



K 



-SF + 2HSF -IH + 3H^- SF 



-F a 



3 a2 
2F + 2HF + gij4>'^''] a 



F 



where i(<j&) = \{tab + ha) and t[ab] = ^{tab - ha)- hab = 
gab + naUb is the projection tensor based on n^. 0, ffab, 
and da are the expansion scalar, shear tensor, and the ac- 
celeration vector based on n^, respectively. The vorticity 
tens or o f the normal - vec tor, LOab natur ally vanishes, see 
eq. (C5). From eq. (C2) using eqs. (^ CI 13 ij) we can 
show 

§ = 3H- K, 

na-w ,,_i„(3)A , fl^(t.)y^(«) ,„2-(t)^(t) 

Oa — OL^a- (C3) 
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Therefore, — k and x can be interpreted as the per- 
turbed expansion scalar and the scalar part of shear of 
the normal-frame, respectively. The trace of the extrin- 
sic curvature is equal to minus of the expansion scalar, a 
and [3 can be seen as perturbations in the lapse function 
and shift vector, respectively. ^''■"^ and c^*^ also cause 
the shear in the perturbed normal hypersurface. 

In order to interprete the velocity related quantities 
we introduce frame-invariant combinations of the four- 
vectors as WA 



(C4) 



Similarly as in eq. (C2) we can introduce the kinematic 
quantities based on iia vector 



'ab = hlh^U(^c-d) — "(a;6) + 0-{aUb), 



<^ab = &ab 
^ab 



6" c;d , 



Qa = Ua-.bu'', 



(C5) 



where hab = gab + UaUb is the projection tensor based on 
Ua- 0, Gab, i^ab, and Oq are the expansion scalar, shear 
tensor, vorticity tensor, and the acceleration vector based 
on Ua, respectively. From eq. (C5) using eqs. (|C^,p 13 1^ 
,0) we can show 



Uq = Uq, Ua 



6 = iH-K+-v^'\ 
a 

^ (3) A 



1/3 



k 



(C6) 



Cabcd — Rabcd — {9acRbd + QbdRac — 9bcRad — 9adRbc) 



R . , 

+7" ygacQbd - 9adgbc) ■ 

6 



(C7) 



It is decomposed into the electric and the magnetic parts 
as 



Eab = Cacbdu'^U , Hab 



,:Vac'fCefbdU%''. (C8) 



Both are symmetric, tracefree and orthogonal to 
Eab = Eba, = = Eabu'', and the same for Hab- 
The nonvanishing electric and magnetic parts of the Weyl 
curvature are: 



E, 



Q/3 



1 



2 



it) + 



A-2K 



Y. 



a/3 ' 



/3)7<5 



which follow from the Riemann curvature tensors and 



eqs. (|5Ul4|). 

In the ADM notation 

500 = <?Oa = A^a, <7a/3 = /la/3, (CIO) 

where is based on /i^^ with the hf^^ the inverse metric; 
only in the rest of this Appendix hap indicates the ADM 
three-space metric. The normal four- vector is tiq = — A^ 
and Ua = 0. The extrinsic curvature is 



K = h"f^K, 



aP, 



(Cll) 



where a colon ' : ' indicates a covariant derivative based on 



hal: 



Pi 



is the connection based on hap- The ADM 



fluid quantities are 

E = UaHbT"-'', Ja 
— uaP < 



-UbTa, Sap = Tap, 



1. 



S = Sap, Sap = Sap — -hapS. 

Compared with the perturbed metric in eq. 

haP = {gal + 2Ca/3 



(C12) 
we have 



N^a{l + A), Na 



-a^Ba, 



r 



(3)c 



Pj- Pi' 

Thus we can show 



Pi 



(C13) 



and similarly for the kinematic quantities based on the Kap — —a 
individual fluid four- vectors U(^i)a- 

Weyl curvature tensor is introduced as 



-gl'^il - A) + B^a\p) + Cap + 2-CaP 



K = 

a 



3-(l-A)+i?"|„+C 
a ' 



= -3H + K = -e. 



R 



(h)a 



Pl& 



R 



(3)a 



-a 



PlS 

I" 



P\Sl 



^p\is 



C. 



i\pi 



'pS 1 ^/37 5' 

= ^ - 4(A + iK)^] , 

E = -To = /i, Ja = aT^ = qa + {^J.+ p)ua, 

S = 3p, S1^ = TT^^)^, (C14) 

Jh)a 



where the intrinsic curvature R 



PlS 



is a Riemann cur- 



vature based on hap; K is the sign of three-space cur- 
vature. Thus, if is proportional to the perturbed three- 
space curvature of the hypersurface normal to Ua- 
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FIG. 1. We present the evolutions of the adiabatic density 
perturbations in the corresponding comoving gauges 
for several components, (i) includes the baryon (dot, black), 
CDM (dot-long dash, blue), photon (long dash, red), mass- 
less neutrino (short-dash, green), and massive neutrino (solid, 
cyan). The two figures are (a) fc/ao = O.lMpc"^ and (b) 
fc/ao — O.OlAfpc^^. Also presented is the ip^ (dot-short dash, 
yellow) where i; is the collective fluid velocity. The parame- 
ters are: h = 0.65, flbo = 0.06, fico = 0.5, fi^o = 5.85 x 10"^ 
Jl^o = 3.99 X 10~^, and fi^^o = 0.44. We consider a flat back- 
ground with vanishing A. The absolute value of the vertical 
scale is arbitrary. 
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FIG. 2. We present the evolutions of in various co- 
moving gauge conditions based on fixing the velocity vari- 
ables of the components, i.e., (pu;;) , and ip^: the baryon 
(dot, black), the CDM v'i'(c) (dot-long dash, blue), the 
photon V'l.i^) (long dash, red), the massless neutrino ^v^^^-^ 
(short dash, green), and the one based on the collective ve- 
locity (dot-short dash, yellow) . We consider three different 
scales: fc/oo = 0.001 (top), 0.1 and 0.2Mpc"^ (bottom). For 
fc/ao = O.OOlMpc"^, the baryon, CDM, and the collective one 
are overlapped (top), and the photon and massless neutrino 
are overlapped (bottom). In order to present the behaviors 
in three scales in one frame, we change the absolute scale of 
the amplitude arbitrarily. The parameters are: f2co = 0.94, 
^^i/mO = 0, and the other parameters are the same as in Fig. 






FIG. 3. The evolution of for different scales: 
fc/ao = 0.0001 (top), 0.001, 0.01, 0.2, and O.SMpc"^ (bot- 
tom). The cases of fc/ao — 0.0001 and O.OOlMpc"^ are almost 
overlapped. The parameters are the same as in Figure |^. 
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FIG. 4. We present the evolution of density perturbation 



in the comoving gauge of the massive neutrino, ^Hy, 



,, for 



different scales: ifc/ao = 0.01 (top), 0.1, 0.2, 0.3 and 0.4Mpc"^ 
(bottom). Figure (a) considers the massive neutrino dom- 
inated model with a parameter f2i/„o = 0.94. Figure (b) 
considers substantial amount of the CDM with parameters 
^^co = 0.5 and f2^„,o = 0.44. In Figure (b) we show the 



evolution of CDM 5, 



(C)«(C) 



(dotted, red) as well. The other 



parameters are the same as in Fig. 
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FIG. 6. We present the power spectra l{l + l)Ci of the 
gravitational wave: the temperature Cf® (top, black), the 
cross correlation Cf^ (green), and the polarizations Cf^ 
(blue) and Cf^ (bottom, red). As the initial condition we 
take the scale-invariant {jlt = 0) spectrum and the solution 
with constant amplitude. We normalize the spectra using 
i{l + 1)C?® = 1 for I — 10. The parameters are the same as 
in Figure n. 
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FIG. 5. We present the power spectra £{£ + l)Ci of the 
scalar-type perturbation: the temperature Cf ® (top, black), 
the polarization Cf^ (middle, green), and the cross correla- 
tion Cf ^ (bottom, blue). We take the adiabatic initial con- 
dition with the scale-invariant (ns — 1) spectrum. Figure 
(a) shows the spectra in logarithmic scale, and (b) shows in 
real scale. We normalize the spectra using £{£ + l)Cf ® — 1 
for £ = 10. The parameters are: flco ~ 0.25, Qao ~ 0.69, 
^I'mO = 0., and the other parameters are the same as in Fig. 
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